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This paper addresses the valuation of the Paris barrier options of [C JY] us- 
ing the Laplace transform approach. The notion of Paris options is extended 
such that their valuation is possible at any point during their lifespan. The 
Laplace transforms of [CJY] are modified when necessary, and their basic 
analytic properties are discussed. 

1. Introduction: This paper is the first part of a report on the valuation of the Paris 
form of barrier options proposed in [CJY] . In their standard form European-style barrier 
options come as puts or calls that are activated or deactivated as soon as their underlying 
hits a prespecified barrier level. The new idea of the Paris barrier options for cushioning 
this abruptness is to introduce a systematic delay for these consequences of hitting the 
barrier to become effective. 

To fix ideas consider the Paris down-and-in call option. Suppose the underlying is moving 
downwards and hits the barrier level. This Paris option is then not immediately activated. 
Instead, a clock starts ticking at the very moment the underlying passes below the barrier 
level. It will tick for D consecutive periods in time. Only then, with this delay of D 
periods the option will be activated if the following two conditions hold. First, one is still 
during the life span of the option. Second, during any of these last D periods of time the 
underlying has been staying below the barrier level. To further illustrate this Paris option, 
suppose D periods before its maturity the option is not yet activated and no clock has 
started ticking. Then the option will not be activated during its life time. It is already 
worthless. 

This paper addresses foundational questions in valuating Paris options. For instance, a 
principal difficulty had been pointed out to the author by Pliska. Suppose the clock is 
already ticking for a Paris option, how should one valuate this option? The necessary 
extensions of [CJY] are discussed in Part I. For valuating Paris options the Laplace trans- 
form approach is adopted. Indeed, Laplace transforms of their prices have been computed 
in [CJY], and the discussion of Part II is based on their results. However, pursuing Pliska's 
suggestions we found it necessary to modify one of these two Laplace transforms. So I am 
very grateful to Pliska for his remark. In this conncetion I also wish to thank Yor for 
making available to me corrections to [CJY] collected by Jeanblanc-Picque. 

This status quo ante set by [CJY] is extended in Part III where a number of insights into 
the analytical structure of the above Laplace transforms are discussed. More precisely, 
these Laplace transforms are constructed using a function \E' whose analytical properties 
are analyzed. In particular the series of this function is derived and also an asymptotic 
expansion for it is proved that is uniform on the right-half plane. These results are 
fundamental for the Laplace inversion of the above Laplace transforms. 

Recall that the Laplace transform generally consists of two stages. As a first step, a suit- 
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ably nice function on the positive real line is transformed into a complex-valued function. 
This Laplace transform exists on a half-plane sufficiently deep within the right complex 
half-plane and defines there a holomorphic, i.e., complex analytic, function. As a second 
step, the Laplace transform has to be inverted to give the desired function on the positive 
real line. The natural methods in this step are complex analytic. Relevant concepts for 
the endeavours ahead are recalled in Part 0. 

To proceed, recall how for valuating barrier options the relative position of strike price 
and barrier play a role. For Paris options this list has to be extended and we propose the 
following notions. If the barrier is equal to or below the strike price we have a standard 
situation, otherwise we have a perverse situation. Also the relative position of today's price 
of the underlying to the barrier level matters. We have the first situation if today's price 
of the underlying is equal to or above the barrier level. If today's price of the underlying 
is below the barrier level, so that the clock for Paris option has started ticking, we have 
the second situation. These situations combine into four possible cases. In each of these 
cases there is a different Laplace transform to be inverted. 

With this program in mind the reader now may wish not only to take a look at this paper 
but also at its companion [SP]. The focus of this last paper is on Laplace inversion. It 
addresses valuating the Paris options in the above first standard case. However, its results 
about Laplace inversion are typical and fundamental for the remaining three cases. 

Acknowlegements: In addition to Professors Pliska and Yor already mentioned in the intro- 
duction, I would like to thank in Mannheim Professor R. Weissauer, Dr. U. Weselmann, and 
Dr. D. Pulea for their interest in and support of this project. 



2. Certain classes of higher transcendental functions: This paragraph recalls 
pertinent facts about the the gamma function and the error function from [L] , [D] . 

Gamma function: The gamma function P is defined on the right complex half-plane 
by the improper integral: 



for any complex number z with positive real part. It interpolates the factorial in the sense 
that r(n +1) = n!, for any non-negative integer n, and more generally satisfies the recur- 
sion rule r{z+l) = zr{z). The functional equation T{z) ■r{l—z) •sin(7rz) = tt in particular 
shows how it continues to a meromorphic function on the whole complex plane with simple 
poles in the non-positive integers as its singularities. It satisfies the duplication formula 



For any complex numbers a, z with positive real part, one has the identity: r{z) = 
Io° e~°'*t^~^ dt. Indeed, for a real it is obtained by a change of variable. Applying the 
identity theorem the resulting identity is then valid on the whole right half-plane. 



Part 



Epitome of methods 




2^'-^T{z)T{z+l/2) = 0FP(2^). 
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Unlike the gamma function itself, its reciprocal is holomorphic on the whole complex plane. 
A repesentation as a contour integral is given by the Hankel formula: 



that is valid for any complex number x with positive real part, and now for any complex 
number a. Herein C is any of the Hankel contours Cg^ji with < £ < i?, or their limits for 
e converging to zero, defined as the following paths of integration: enter from minus infinity 



on a parallel to the real axis through —ie until circle in the orgin with radius R is first hit. 
Pass counterclockwise on this circle until the parallel to the real axis through +is is hit in 
the left complex half-plane. Leave on this parallel to minus infinity. 

Error function: The complementary error function Erfc is the analytic function on the 
complex plane given by: 



for any complex number z, where the path of integration in the last integral starts at z 
and moves parallel to the real axis with the real part of its elements going to infinity. 

3. Laplace transform: This paragraph collects pertinent facts about the Laplace 
transform from [D]. 

The class of functions considered in the sequel is that of functions of exponential type, i.e., 
of continuous, real-valued functions / on the non-negative real line such that there is a real 
number a for which exp {at)f{t) is bounded for any t > 0. The Laplace transform is the 
linear operator L that associates to any function / of exponential type the complex-valued 
function L(/), analytic on a suitable complex half-plane, given by: 



for any complex number z with sufficiently big real part. The convergence theorem is 
that the Laplace transform L(/) converges on the whole half plane {2|Re(^) > Re (2:0)} 
if L(/)(2o) exists. So there is a unique —00 < cr = cr/ < 00, the abscissa of convergence 
of L(/), such that L{f){z) converges, for any z with Ke{z) > a and does not exist for z 
with Re [z) < a. To fix ideas, consider the function given by f{t) = exp (at) with a any 
positive real number so that formally 'L{f){z) = Jq° exp {—{z — a)t) dt. With the absolute 





Figure A The Hankel contour C^^ji 
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value of cxp being equal to exp(Rc(u;)), for any complex number w, this improper 
integral exists iff the real part of z — a is positive, i.e., iff the real part of z is bigger than 
a. In this case, Ij{f){z) equals {z — a)~^ and is analytic on the half plane {^|Re (z) > a}. 

That L is an injection is a consequence of ultimately the WeierstraB approximation theo- 
rem. Its inverse, the inverse Laplace transform L~^, is expressed as a contour integral by 
the complex inversion formula of Riemann. For any function H analytic on half-planes 
{Re (z) >xo} with xq any sufficiently big positive real number, it asserts: 

L-\H){t) = — / e''-H{z)dz, 

for any positive real number t if H satisfies a growth condition at infinity such that the 
above integral exists. In particular (L~^ o L)(/) = / if / is of exponential type. 

The convolution of two functions / and g of exponential type is the function of exponential 
type given by f*g{t) = Jq f{u)g{t — u) du for any t > 0. With du being Lebesgue measure, 
notice f * g = g * f ■ The product theorem for the Laplace transform is that L transforms 
convolutions into products: 

Uf*g)iz) = -Lif)iz)-Mg)iz), 

for any z with Re{z) sufficiently big. Equivalently, / * g{t) = L~^(L(/) ■ 'L{g)){t). The 
proof is a formal Fubini calculation as soon as one extends functions of exponential type 
by zero on the non-positive real line, whence / * g{t) equals f{u)g{t — u) du. 

There are a number of shifting principles for the Laplace transform. First, weighting any 
function / of exponential type with the function sending any t > to exp (at), where a is 
any real number, introduces a shift by minus a in the Laplace transform: 

L(e«V)(^) = L(/)(z-a). 

The effect of weighting the Laplace transform of / by exp {az) depends on the sign of the 
real number a. More precisely, let a be any non-negative real number. Using the above 
complex inversion formula, one has: 

L-i(e-L(/)(^))(t) = /(t + a), 

for any t > 0. Weighting with exp {—az) gives the following shifting theorem: 

L-i(e— L(/)(^))(0 = l(,,^)(t)/(t-a), 

for any positive real number t. Herein \x is the characteristic function of any set X. This 
is proved calculating the Laplace transform of the right hand side. 
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The integration and differentiation principles for the Laplace transform describe the be- 
haviour of the Laplace transform upon multiplication by integral powers of the complex 
variable z. Using partial integration in the defining integral of the Laplace transform: 

L(/('^))(2) = 2^ ■ L(/)(2) - /(o)(0)z"-i /(^-i)(0)20, for any non-negative integer n 

and any complex number z in the half-plane of convergence for L(/). In particular, 

L(/(-))(^) = ^--L(/)(^), 

if / is such that for any non-negative integer k less than n its k-ih. derivative f^^^ vanishes 
at zero. This last formula holds in general if one considers n-fold integration as differ- 
entiation of the negative order — n. More precisely, the Laplace inverse of z~'^ ■ Ij{f){z), 
for any non-negative integer n, is given as the n-fold iterated integral In{f) of /• This is 
recursively defined by /o(/) = f i and 

Im+l{fm= f Im{f){u)du, 

Jo 

for any positive real number t, and any non-negative integer m. 

Part I Paris options 

4. Black— Scholes framework for valuating contingent claims: For our analysis 
we place ourselves in the Black-Scholes framework and use the risk-neutral approach to 
valuating contingent claims. In this set-up one has two securities. First there is a riskless 
security, a bond, that has the continuously compounding positive interest rate r. Then 
there is a risky security whose price process 5" is modelled as follows. Start with a complete 
probability space equipped with the standard filtration of a standard Brownian motion on 
it that has the time set [0, cxo). On this filtered space one has the risk neutral measure Q, 
a probability measure equivalent to the given one, and a standard Q-Brownian motion B 
such that S is the strong solution of the following stochastic differential equation: 

dSt = ir-6)- Sfdt + a- Sf dBt, t e [0, oo). 

Herein the positive constant a is the volatility of S. The constant d depends on the security 
modelled. For instance, it is zero if is a non-dividend paying stock. 

5. Paris options: This paragraph dicusses the Paris barrier options of [CJY, §2]. 
Their new idea for cushioning the impact of the underlying hitting the barrier is as fol- 
lows. They require the underlying S to spend a minimum time D > above or below their 
prespecified barrier L > before the option is knocked in or knocked out. 

The Paris down-and-in call to be considered in the sequel is the following European-style 
contingent claim on S written at time to and with time to maturity T. Its payout at T is 
that of a call on S with exercise price K: 

(St - K)+ = maxj^T - i^, 0} , 

if St is less than L during a connected subperiod of length at least D of the monitoring 
period [to, T]. Equivalently, there is a point in time a such that the interval / = {a,a + D) 
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is contained in the lifespan [to)2^] of the option and St < L for any t in /. Otherwise the 
call expires worthless. 

In the sequel, fix any time t during [to)?"] at which the Paris down-and-in call has not 
yet been knocked in. The time-t value Cd,i of this call is then a function of K, t = T — t, 
So, r, S, L and D. In the sequel, mostly any of these arguments is suppressed from the 
notation. The valuation problem depends on an excursion time Hl ^ — H£ j-, ^ defined in 
the sequel such that the Paris down-and-in call is knocked in before its maturity T if and 
only if HL,t < T. Using the arbitrage pricing principle, the time-t price Cd,i of the Paris 
down-and-in call is then given by the conditional Q-expectation: 



0(5't) • l{H^,t<T}], 



with (f){ST) being equal to the plain vanilla call payout {St — K)^ , and denoting the 
characteristic function of the event Hl i < ^ by 1{Hl t<T}- 

In defining Hl i at time t two constellations are to be distinguished. Namely, either to- 
day's security price St is equal to or above the barrier L, or it is below L. In the first case 
with St> L, the Paris down-and-in call is knocked in before its maturity T iff S' remains 
smaller than L for all points in time of an interval of length at least D contained in [t, T]. 
Thus let Hi t in this case denote the first time s greater than or equal to t-\-D such that 
Su < L for any time u in the interval {s — D,s). 

In the second case where St < L, the security price S has already been staying below L 
for some connected period time during the lifetime of the option. The down-and-in call is 
then knocked in before its time to maturity T if the following happens. The security price 
S continues to stay below L also during the period of time from today until the point in 
time t + d with d smaller than D, i.e., one has Su < L for all u in {t,t + d). Thus define 
H^ t = t + d in this case. If the level L is hit by S before time t + d, however, the clock 
for the minimum lenght D is restarted. So define i^i,,t in this case to be the first time s 
greater than or equal to t+D such that Su < L for any u in the subinterval {s — D,s) of 
the positive real line. 

The discussion of Paris barrier options would not be complete without having referred to 
[CJY, §2] for details on the whole family of Paris barrier options. Calls and puts with ei- 
ther of the following barrier types: down-and-out, down-and-in, up-and-out, up-and-in. 
Valuating these eight types of options is reduced in a standard way to the case of the above 
Paris down-and-in call. Indeed, put-call parities are given in [CJY, §5], and [CJY, §4.2] 
reduces the valuation of the out-calls to that of the in-calls, mutatis mutandis. This paper 
thus concentrates on valuating the Paris down-and-in call, in the sequel also referred to 
as the Paris option for simplicity. 



Part II Laplace transforms 



6. Laplace transforms of the value of the Paris option: This paragraph discusses 
the Laplace transforms of the densities for the the Paris option of paragraph five. 
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The basic valuation identity: This notion of density is made precise by by the follow- 
ing basic valuation identity for the time-t price Cd,i of the Paris option: 

Cd,i^e~^'^^>'" / e^^ ■{Sfe''^-K)-h^{r,x)dx, 

with the Paris option density hb given by: 



hb{u,y) = E 



oo 



oo 



. {x-yf 



dii*{x) , 



for any real numbers u > and y, and where P{St) — log{K / St) , where = HL,t — t, 
and with fx* = lJib,D,- the measure for Brownian motion at time H^. By construction, 
hb{u, ) is zero for u less than D, respectively d depending on today's situation. 

Laplace transforms of the densities: With the proof of the basic valuation identity 
postponed to the end of this paragraph, the Laplace transforms of with respect to time 
are discussed next. For any fixed real number y they are recalled to be defined by: 



/•oo 

L(/ib( ,y)){z) = / e-''hb{t,y)dt, 
Jo 



for any complex number z with sufficiently big real part. Moreover define the function ^ by: 

/•oo ^ 

^i^) — X ■ e ^ dx , 
Jo 

for any complex number z, and choose on the complex plane with the non-positive real 
axis deleted the square root defined using the principal branch of the logarithm. The 
computations of [CJY, §§5, 8] can then be summarized as follows. 

Proposition A: For any real number y, the Laplace transform of hb{ ,y) is a holomor- 
phic function on the right complex half-plane. 

Proposition B: Suppose b = a~^log{L/St) is non-positive. For any fixed real number 
y, the Laplace transform of hb{ ,y) is given by: 



eVD r°° _i 



\h—x — y\\/2z 

for any complex number z with positive real part. 



If b = a~^log{L/St) is positive, i.e., when today's security price St is below the barrier 
L, we have modified the Paris option of [CJY]. Today's situation is then such that the 
price of the security has already been staying below the barrier L for a certain connected 
period of time during the lifetime of the Paris options. For the Paris option to be knocked 
in the price thus has to stay below the barrier L only for another connected period time 
of a length d < D from today, i.e., S^ < L for all u in [t, t+d). In this way, Paris options 
now can be valuated at any point of their monitoring periods. The desirability of this 
modification has been pointed out to me by Pliska. I regard the following modification of 
the valuation results of [CJY] as a direct result of his suggestions. 
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Proposition C: Suppose b — a~^log{L/St) is negative. For any real number y, the 
Laplace transform of hb{ , y) is given on the right half-plane as the following four-term 
sum of Laplace transforms: 

Hhi ,y)) =Erfc(-|=) .L(;jt,i( +L(fti,,2( ,!/)) 



Herein hf^^k are the functions on the positive real hne times the real line defined by: 

1 r / T fa-'')' 

h r \ 1 / \ 1 J Vu-d 



2y/u-d+D \y/2{u-d){u-d+D) J 

hhA{u,y) = l(d,oo)(«)^== ■ fb,uiy) where fb,uiy) = e ^ - e , 

for any positive real number u and any real number y, while the remaining two functions 
are defined using their Laplace transforms as follows: 

1 e~^^ ( f ^/2~ \ 



V27Td Jo \pl^-^{^2Dz) Via 

for any positive real number y and any complex number z with positive real part. Herein 
\X}) denotes the law of the first passage time to the level 6, given by 

b / --^ 

IJb{dw) ^ —=w~^^'^e 2"-, 



/2tt 

on the positive real line. 

Proof of the basic valuation identity: The proof of the basic valuation identity is 
by reduction to that of [CJY, §§4, 5]. As a first step, the valuation problem for the Paris 
option of paragraph two is transcribed in in terms of a restarted-at-time-t normalized 
Brownian motion with drift W* as follows. With B* the restarted-at-time-t Brownian 
motion B*{u) = B{t+u)—B{t), consider the Brownian motion with drift: 

W: = i log(^) = wu+Bl , I. e [t, oo), 
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whose drift coefficient is zu = ■ {r — 5— cr^/2). Recalling h = ■ log{L/St), one has 
St+u < iff W*{u) < b hj construction. Thus define the excursion time = H^'^ q for 
W* as follows. If today's value VF*(0) of W* is greater than or equal to b then is the 
first time s greater than or equal to D such that W* {u) < b for any time u in the interval 
(s — D^s) of the positive real line. Suppose today's value W*{0) of W* is less than b. Set 
= dif W*{u) < b for any time u in the interval [0, d) remaining for W* to have stayed 
a connected period of time of length D below L. Otherwise, is the first time s greater 
than or equal to D such that W*{u) < b for any time u in the subinterval {s — D,s) of 
the positive real line. One has = HL^f-t, whence t <T iS <r = T—t. Using 
Strong Markov, the valuation problem for the Paris option thus transcribes as follows: 

+ 



d,i 



St 



-K) 



As a next step, using the particular case [KS, p.l96f] of the Cameron-Martin-Girsanov 
theorem, change measure from Q to Q* such that the drift term in W* is killed and W* 
becomes a (5*~Brownian motion. Change measure in the expectation expressing the price 
Cd,i and iterate the Q*-expectation thus obtained to get: 



E* 



MH;<r}E* 



(St 



■K) 



+ 



Using as F-stopping time apply the strong Markov property of Brownian motion in 
the conditional expectation. The resolvent thus obtained can be further simplified with the 
random variables ELI ^'^'^ being independent. With //* the measure for 

the basic valuation identity follows. 

7. An intermediate identity for the proofs of the Laplace transforms: Adopt 
the notation introduced in the previous paragraph six in the basic valuation identity for 
Paris options. Then one has the following intermediate identity: 



/ E* 




e 


In 





■-\x—y\V2z 



Lemma: For any real number y, one has: 



for any complex number z in the right complex half-plane. 



II* {dx) , 



For proving the Lemma interchange the Laplace transform with the two exponential inte- 
grals defining the Paris option density y) in paragraph six to get: 



L(/i6( .y)){z) 



E* 



R 



1 



fi*{dx) . 



In this last integral successively change variables w = u—H^ and separate the expectation 
from the Laplace transform to get: 



/ E* 













1 



{x-yf 
1u 



.^/2 



TTU 



(z) IJL*{dx). 



Using [Sch, §3] for the Laplace transform completes the proof of the Lemma. 
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8. Computing the Laplace transform in the case b equal to zero: The compu- 
tation of the Laplace transforms of paragraph six for is by reduction to the case b equal 
to zero where one has to show the key relation: 



E* 



^{zVD) = ^{0) = 1, 



for any complex number z with positive real part. This paragraph reviews the key ideas 
of its proof in [C JY, §8] . The computations in this case b equal to zero are based on the 
Azema martingale and properties of the Brownian meander. The Azema martingale is the 
martingale on [0, oo) for the slow Brownian filtration F"*" given for any t > by: 

Ht = isgnW^)^/t-gt. 

Herein is defined as the supremum over all real numbers s < t such that W*{s) — 0. 
The process given by: 



mt{u) 



IW'^igt + uit-gt))], «e [0,1] 



is a Brownian meander and is independent of the a-subfield F~^{gt) of the slow Brownian 
filtration. Its law is independent of t and given by a;exp {—x^/2) ■ l(o,oo) dx. The tautology 
W*{t) = mt{l) ■ jit thus implies the following identity: 



zW* 



— e ^ 



for any positive real number z. Fixing any such z, the crucial point is to exhibit the product 
iI){—z-ii{Hq)) exp {—z^Hq/2) as a martingale. This is achieved by a boundedness argument 
showing the argument of the above conditional expectation to be uniformly integrable for 
any time t up to Hq . Using another form of the optional stopping theorem it follows that 
the (5*-expectation of the martingale stopped at Hq equals the Q* -expectation of it at 
time zero, whence ^(0) which equals one. With the random variables Hq and W*{Hq) 
independent, the key relation results. This is an identity between functions holomorphic 
in particular on the right half-plane. Using the identity theorem, it thus remains valid for 
any complex number z with positive real part. 



9. Computing the Laplace transform for b non— positive: Using §7 Lemma the 
proof of §6 Proposition B reduces to compute the expectation of exp {—zH^), for any com- 
plex number z with positive real part, and the density //*. Following [CJY, §8.3.3], this 
is by reduction to the case b equal to zero of the previous paragraph. Indeed, decompose 
as follows: 

Herein Tb is the first passage time of W* to the level 6, and Hq* is defined as follows. It 
is the smallest non-negative point in time s at which the restarted-at-time-Tb Brownian 
motion W**{u) := W*{Tb + u)-W*{Tb) = W*{Tb + u) - bis zero for the first time after 
having been less than zero for a connected period of time of length at least D. 

To compute the expectation required, notice that T5 and Hq* are independent random 
variables. The conditional expectation at time T5 of exp {—{w'^/2)H^) thus is the product 
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of exp {-{w"^ /2)Th) and the expectation of cxp {-{w"^ /2)Hq*). With the law of Hq* equiv- 
alent to that of i^Q, this last expectation is given by the key relation of paragraph eight. 
So it is deterministic in particular. Take expectations to time zero of this product. The 
expectation of exp ( — (t(;^/2)T6) that remains to be computed is the Laplace transform of 
the law of % at Using [Sch, §3] it is equal to exp(-|6|(2 • {w'^/2)Y^'^) and thus to 

exp {-\-hw). With w = (22;)-'^/^ one gets: 



■m 



^{V2D^ ) 



E* 



^{V2D^ ) 



To determine /x* notice the tautology W*{H*) = W**{H**) + b. Thus W*{H*) < x iS 
W**{Hq*) < x — h. From the discussion of the case h equal to zero, the law of W**{Hq*) 
is that of the negative of mi . As a consequence: 



Q*{W*{Hl) e dx) = l(_oo,6](a^) ■{b-x)-e 



2D 



dx 

'd 



Substitute into §7 Lemma and change variables w — b — x to obtain §6 Proposition B for 
any complex number z with sufficiently big positive real part. Extend to the whole right 
half-plane using the consequence of the key relation of paragraph eight that ^{{2Dz)^/'^) 
has no zeroes there. This completes the proof of §6 Proposition B and the 6-non-positive 
part of §6 Proposition A. 

10. Computing the Laplace transform for b positive: The proof of Proposition C 
of paragraph six modifies the argument of [CJY, §8.3.3]. Using §7 Lemma, the first task is 
to compute the density n* and the (5*-expectation of exp {—zH^) for any complex number 
z with positive real part. The two cases where the level b is hit before or after time d are 
to be distinguished. Thus decompose the underlying probablity space 0, into the set A 
on which the first passage time Tb of W* to the level b is less than or equal to d and its 
complement $7 \ on which Tf, is bigger than d. This induces the decomposition: 



E* 



= E* 



U • e 



+ E* 



of the Q* -expectation of the random variable ex.-p{—zHl). By construction — d on 
the complement of A. So the second above summand is e'x.\){—zd)Q*{Tb > d). On the set 
A the level b is reached before the critical time d and the clock for the excursion is reset 
to zero at the first passage time T5 to the level b. Accordingly one has the decomposition: 

with Hq* defined as follows. It is the smallest s > at which the restarted-at-time-T;, 
Brownian motion W**{u) := W* {Th+u)—W* (Tb) = W*{Tb+u)—b is zero for the first time 
after having been less than zero for a connected period of time of lenght at least D. 

With this random variable Hq* independent from T^, the conditional expectation of the 
random variable 1a ■ exp {—zH^) at time thus equals exp {—zTb) times the expectation 
of ex.p{—zHQ*). The key relation of paragraph eight now applies to Hq* and identifies 
this last expectation as the reciprocal of ^((2zD)^/'^). The expectation of 1a ■ exp (—zTb) 
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on the other hand, is given by integrating exp (—zw) from zero to d against the density 
fXb of Tft. Summarizing, it so foUows for any complex number z with positive real part: 



For determining //* decompose the random variable W*{H^) with respect to A: 

W*{Hl) = {W* * iH**)+b) ■ 1a + W*id) ■ ln\A ■ 

Using the independence of T;, and Hq* , the law of the first summand is obtained as the 
convolution of the laws of Tb and W**{Hq*), whence 



Ulidx) = Q*{{W*%H**)+b) e dx;n < d) 

dx 

'd 



Q*(Tb<d)-l(_oo,6i(x)-(6-x)-e ^ 



For the law of the second summand notice T5 > iff W* (t) < b for all t < d, ot equivalently, 
max{W* {t)\t < d} is smaller than b. Using [H, p. 9] one so has: 

li^idx) = Q*{W*{d) e dx;Tb > d) 

= e'^ .^^-e~^^dx■ 



V2^Td Vm ' 

With the expectation factor in §7 Lemma seen above to be deterministic and independent 
of the variable x, the Laplace transform of hb{ , y) at z is given as the product: 

L{hb{,y)){z)^E*[e-^'''^] [ ' ' ^ {i^l + i.*,){dx) . 

Jn \J 2.Z 

One is so reduced to compute the two integrals of the second factor of this product. For 
the first of these, on substituting for //* and changing variables w = b—x, one gets: 



L 



lii{dx) = — / X- 7= -e dx . 



For the second of these one analogously obtains: 



L 



-\x-y\V2z \^-y\^2z _ix-2bf , 

, = — ii*(dx) = ^= / = — (e 23 _e ]dx. 

At this point it remains to identify the two Laplace inverses of §6 Proposition C to complete 
its proof. This is based on the result: 
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valid for any positive real numbers a and u, and discussed in [Sch, §3]. Applying Fubini's 
theorem to get the inversion integral as inner integral, the improper integral factor of the 
above integral becomes: 



^0 



oo ^2 / ^-ib-x-yW^z 



xe 



^ • L-M e-"^- = ){u) dx 



Successively apply the shifting theorem for the Laplace transform of paragraph three to 
take care of the factor exp {—zd) and the above Laplace inversion formula with a equal to 
\b—x—y\ to arrive at: 

1 1 r°° ( I (b-x-yf \ 

D ' ^2T,{u-d) Jo 

This integral is seen to be the value of the function hb,3 at u and y. Analogously, the Laplace 
inverse at any positive real number u of the integral factor of the above iJ,2 integral is: 

1 r -i^niiL / _xL {x-2bf ^ 

e 2(u-^) fe 2d _e 2d ]dx 



and is seen to be equal to the value of the function hb,4 at u and y. This completes the 
proof of §6 Proposition C. 



Part III Analytic properties of the function ^' 

11. Analytic properties of This paragraph collects and discusses pertinent ana- 
lytical properties of the function "if for any complex number w recalled to be given by: 

POO ^.2 ^ 

^{w) = X ■ e ~^ dx . 
Jo 

Developing the linear exponential factor of the integrand in its series and integrating the 
resulting series term by term, gives the following series expansion: 

^(w)=y — -r — — --u;^. 

^ Z^n=Q n\ V 2 / 

This series is absolutely convergent for any complex number w, and convergence is uni- 
form on compact sets. Rearrange it in its even and odd order terms. Replacing w by 
its negative leaves unchanged the even part and produces a minus sign in the odd part. 
Using the duplication identity for the gamma function and redeveloping the square order 
exponential series that results as a factor, then gives the following key identity: 

= *(-'«;) -I- V2tt • w • e 2 
connecting the values of on the right half-plane with those on the left half-plane. 
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12. A uniform asymptotic expansion: The aim of this paragraph is to prove for 
the function \1/ given by = Jq xexp + wx) dx, for any complex number w, 

the following uniform asymptotic expansion. 

Lemma: For any complex number w with positive real part, one has the expansion: 

^ fe=i 

for any positive integer N, whose remainder term Rn+i satisfies the estimate: 

(2iV+l)! 1 



\Rn+i{w)\ <2 



Corollary: For any complex number z not a non-positive real number, one has: 

*(-^ ) = - +i?2(^ ) 
Z 

with i?2(-z^^^) being big Oh in \z\~'^ for the absolute value of z going to infinity. 

The Corollary is an immediate consequence of the Lemma. Recall that the principal branch 
of the logarithm has been chosen. The proof of the Lemma is in two steps. As a first step, 
establish the asymptotic expansion with the remainder term satisfiying the estimate: 

IP f M ^ Cm 1__ , „ (2iV+l)! 

if I arg(w)| < n/ 2 — 6yj lor 6yj > 0. For this develop the square exponential factor of the 
integrand for '^{—w) into its Taylor series up to the {N—l)-st term and integrate the 
resulting sum term by term. Using the duplication formula for the gamma function the 
coefficients of the asymptotic expansion follow on shifting the summation index by one. 
Proceeding analogously, the absolute value of the remainder term can be majorized hy Cn 
times the reciprocal of the 2(A/'+l)-st power of the real part of w. With w in the right 
half-plane. Re (w) equals \w\ ■ cos{9) for an angle 9 of absolute value less than or equal to 
t:/2—6w as above. Thus cos(^) is bigger than the positive constant cos(7r/2— ^^y) completing 
the proof of the first step. 

As a second step, extend the above asymptotic expansion by showing the above estimate 
to be valid also for any complex numbers w such that 7r/4 < | arg{w)\ < Stt/A — Sw This 
is based on the following integral representation: 

Jo 

for any complex number w. Herein 9 is any fixed any angle of absolute value less than n/4 
and the integral is over the path sending any positive real number RtoR- exp {i9). This 
representation is proved using Cauchy's theorem. Indeed, the absolute value of the above 
integrand is majorized by |^| • exp (—Re (C^)/2 — Re (w^)). It thus converges uniformly to 
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zero with the radius R going to infinity on that part of the circle in the origin of radius R 
which is parametrized by angles of absolute value less than 7r/4. 

Now fix any positive real number s smaller than Sw If arg(w) is negative, choose the angle 
9 = +n/4—e in the above integral representation. Otherwise let 9 = —n/4+e. In the sequel 
the case axg{w) positive is considered. In the above integral representation then develop 
the square exponential factor of the integrand into its Taylor series up to the {N — l)-st 
order. In the identity: 

Jo Jo 

write wexp{i9) = exp (z(6' + arg(t(;))). Notice that the angle 9 + airg{w)) is between 
9+n/4: — e and 9 + ?)n / A — dyj — 7r/2 — {d^ — e). Thus it is smaller than 7r/2 in particular. 
As a consequence, the real part of = wexp{i9) is positive. Recalling that r(a + l) is 
equal to the integral Jq exp {—r]t)t°' dt, the above integral is thus seen to be equal 
to . r(2(/c + l)). The coefficients of the asymptotic expansion follow. Proceeding 

analogously, the absolute value of the remainder term Rn^i{w) is majorized by Cjv times 
the reciprocal of the 2(A^+l)-st power of \w\ ■cos(6'+arg(2;)). Recalling the above bounds for 
its argument, the cosine herein is positive and bounded from below by cos(7r/2 — (5^^— e)), 
for any positive e less than §yj. Thus it is bounded from below by cos(7r/2— (5ty), completing 
the proof of the second step. 

To complete the proof of the uniform asymptotic expansion, apply the estimate of the first 
step with any Sg strictly between 7r/4 and 7r/2 if |arg(t(;)| is less than 7r/4. Otherwise, 
apply the estimate of the second step with 5yj strictly between 7r/4 and 7r/2. In both cases, 
cos(7r/2— is minorized by cos(7r/4) = 2~^/^. Cancelling the powers of two one gets the 
factor two. This completes the proof of the uniform asymptotic expansion. 



References: 

[CJY] M. Chesney, M. Jeanblanc-Picque, M. Yor, Brownian excursions and Parisian barrier 

options, Adv. Appl. Probability 29(1997), 165-184 
[D] G. Doetsch, Handbuch der Laplace Transformation I, Birkhauser Verlag 1971 
[H] J.M. Harrison, Brownian motion and stochastic flow systems, Krieger reprint 1990 
[KS] I. Karatzas, S.E. Shreve, Brownian motion and stochastic calculus, 2nd ed., GTM 113, 

Springer 1991 

L] N.N. Lebedev, Special functions and their applications, Dover Publications 1972 
Sch] M. Schroder, Fact and fantasy in recent applications of transform methods to the valua- 
tion of exotic options, Fakultat fiir Mathematik und Informatik, Universitat Mannheim, 
September 1999 

[SP] M. Schroder, On the valuation of Paris options: the first standard case, Fakultat fiir 
Mathematik und Informatik, Universitat Mannheim, November 1999 



On the valuation of Paris barrier options: 
foundational results 

Michael Schroder 

(Mannheim) 



November 1999 



